x=-3:0.05:13;
y=3.0308*x.^2-3.0550*x+0.2468;
X=[-3 -1 1 13];Y=[15 5 1 5];
plot(x,y,'--r',X,Y,'*g')